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Abstract 

We calculate the expectation values of the stress-energy bitensor defined at two different space- 
time points X, x' of a massless, minimally coupled scalar field with respect to a quantum state at 
finite temperature T in a flat A^-dimensional spacetime by means of the generalized zeta-function 
method. These correlators, also known as the noise kernels, give the fluctuations of energy and 
momentum density of a quantum field which are essential for the investigation of the physical 
effects of negative energy density in certain spacetimes or quantum states. They also act as the 
sources of the Einstein-Langevin equations in stochastic gravity which one can solve for the dynam- 
ics of metric fluctuations as in spacetime foams. In terms of constitutions these correlators are one 
rung above (in the sense of the correlation - BBGKY or Schwinger-Dyson - hierarchies) the mean 
(vacuum and thermal expectation) values of the stress-energy tensor which drive the semiclassical 
Einstein equation in semiclassical gravity. The low and the high temperature expansions of these 
correlators are also given here: At low temperatures, the leading order temperature dependence 
goes like while at high temperatures they have a dependence with the subleading terms 
exponentially suppressed by e~'^. We also discuss the singular behaviors of the correlators in the 
x' — )• x coincident limit as was done before for massless conformal quantum fields. 

PACS numbers: 



* Email: lhtcho@mail.tku.edu.twl 
tEmail: lblhu@umd.edul 



2 



I. INTRODUCTION 



In this paper we present a calculation of the stress tensor correlators, also known as 
the noise kernels, of a massless, minimally coupled scalar field at a finite temperature T 
in a flat A^-dimensional spacetime by means of the generalized zeta-function method. The 
stress tensor correlator is the expectation value of the stress-energy bitensor defined at two 
separate spacetime points. The zeta-function method was first introduced by Dowker and 
Critchley [l| and Hawking [2i] (see also jsj) and successfully applied to the regularization 
of ultraviolet divergences in stress tensors of quantum fields in curved spacetimes with 



Euclidean sections, e^., static b 
by Phillips and Hu |^ (see also 



ack holes and (anti-)de Sitter space. It was generalized 
5I) for the calculation of stress energy correlators. In an 



earlier paper jol we mentioned three classes of physical problems of current interest which 
necessitate the knowledge of such quantities and motivated us to undertake this task for 
AdS spaces. Here we also mention three classes of problems as motivation for our present 
endeavor of calculating the stress energy correlators for finite temperature quantum fields 
in curved spacetimes. Common to both endeavors foremost is 

n 

A) the semiclassical stochastic gravity program [7], a theoretical framework established in 
the 90s as a natural extension of semiclassical gravity (SCG) theory of the 80s jsl, for includ- 
ing the effects of fluctuations in the quantum matter field through the Einstein-Langevin 
equation (ELEq) which govern the behavior of the induced metric fluctuations [9]. While 
SCG goes beyond quantum field theory in curved spacetime (QFTCST) of the 70s 1^ 
(viewed as the test-field approximation of SCG on a fixed background spacetime), in that the 
backreaction of the quantum matter field on the dynamics of the spacetime is incorporated 
through the expectation value of the stress-energy tensor as the source of the semiclassical 
Einstein (SCEq) equation, stochastic gravity goes beyond SCG in that it includes also the 
backreaction of the fluctuations of the stress-energy tensor, measured by the expectation val- 
ues of the stress-energy bitensor, also known as the noise kernel, which govern the behavior 
of the induced metric fluctuations. 

B) Black Hole fluctuations and backreaction: Owing to the existence of an Euclidean 
section in static (e.g., Schwarzschild) black hole spacetimes quantum field effects can be 
obtained via thermal field methods. E.g., Hawking effect can be derived with the use of 
the Euclidean Green function jll,] whose periodicity corresponds to the inverse Hawking 
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temperature. Fluctuations and backreaction of the stress-energy tensor near the black hole 



event horizon is an important issue 



12| . the investigation of this problem for Schwarzschild 



black holes requires a knowledge of the thermal stress tensor correlator, which is the expecta- 
tion value of the stress energy bitensor with respect to the Hartle-Hawking state. Thermal 
stress tensor correlator calculation in hot flat space is the logical first step towards this 
goal, as was the motivation in the earlier work of Phillips and Hu [3]. In that paper 
the authors use the Gaussian approximation [l4| for the Green function for such quantum 
fields to evaluate the noise kernel in two optical metrics: hot fiat space and the optical- 
Schwarzschild spacetime. The optical metric for an ultrastatic spacetime has the product 
form ds^ = gabdx°'dx^ = dr"^ + Qijdx^dx^ . In the Euclidean sector we can allow the imaginary 
r time to possess a periodicity of 27r/fi; = /3 = 1/T , with T the temperature, to connect with 
thermal field theories. (For a black hole, k is the surface gravity.) For massless conformally 
coupled quantum fields in hot fiat space, the Gaussian Green function is exact. For optical 
Schwarzschild, the Gaussian Green function is known to be a fairly good approximation 
for calculating the stress tensor which involves second covariant derivatives of the Green 
'unction. The noise kernels involve up to four covariant derivatives of the Green function 
isl . The Schwarzschild metric is obtained from the optical Schwarzschild by a conformal 
transformation. 

This earlier work is followed by a recent paper 16|] wherein the authors, instead of seeking 



the coincident limits of the noise kernel, calculated the point-separated expression (in both 
timelike and spacelike directions) necessary to solve the equations of stochastic semiclassical 
gravity, by the same Gaussian approximation for the Wightman function for conformally 
invariant fields. These authors have computed all components of the noise kernel exactly for 
hot fiat space and several components for the Schwarzschild spacetime. They showed that the 
noise kernel for the conformally invariant field has a simple scaling behavior under conformal 
transformations which enables them to obtain the results for Schwarzschild spacetime from 
that for the optical Schwarzschild metric. 

C) Negative energy density and fluctuations Apart from applications to early universe and 
black hole physics, closer to home, the stress energy correlators are directly relevant to the 
fluctuations of vacuum energy density, especially the existence of negative energy density as 
in flat space with boundaries (such as the Casimir energy) and for nonclassical states (such 
as two- mode squeezed states). This issue was raised by Kuo and Ford [l7|, furthered by 
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Wu and Ford [18|] and pursued by Phillips and Hu j^, ll9|] , where the fluctuations in energy 
density are shown to be comparable to its mean value for several classes of spacetimes (e.g., 
Einstein Universe) and states (e.g., Casimir). Many related problems such as the quantum 



interest principle 



20| bear on foundational issues of quantum fleld theory and spacetime 



structure. With generalization to flnite temperature one can ask how these quantum vacuum 
behaviors are altered by thermal fluctuations. Results for hot flat space were obtained in 



13| for conformally invariant quantum flelds. Our present results for massless minimally 



coupled flelds in arbitrary dimensions complement and further these earlier studies. 

The organization of this paper is as follows: In the next section we introduce the gen- 
eralized zeta-function method, in Sec III we use this method to calculate the stress energy 
tensor correlator for a flnite temperature quantum fleld in flat N-dim space. In Sec. IV we 
give the expressions for high and low temperature expansions showing which components 
in the bitensor are dominant and how temperature and dimensionality alter this behavior. 
We also discuss the singular behaviors of the correlators in the x' — )■ x coincident limit. We 
conclude with some remarks in Sec. V. The detailed expressions are collected in Appendices 
A, B and C. Our convention is the same as in our earlier papers 6|. 



II. GENERALIZED ZETA-FUNCTION METHOD 



In this section we flrst introduce the generalized zeta-function method of Phillips and 



^ (see also jsl) based on the original methods of Dowker and Critchley |l| and Hawking 



Hu 

a. 



This method has recently been used to consider the case of the stress-energy correlations in 
AdS spaces |6|. Here we apply it to calculate such correlations in hot flat spaces. 

To be general we start by considering a massive m scalar fleld (p coupled to an A^- 
dimensional Euclideanized space (with contravariant metric g^^{x), determinant g and scalar 
curvature R) with coupling constant ^ described by the action 



(1) 



where H is the quadratic operator 



H 



(2) 



and R is the scalar curvature. The effective action deflned hj W = InZ is related to the 
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generating functional Z by 



Z 



Vcpe 



(3) 



The expectation value of the stress-energy tensor can be obtained by taking the functional 
derivative of the effective action 



5W 



(4) 



This formal expression is divergent at coincident limit and some procedure of regularization 
need be implemented. Here we adopt the (^-function regularization scheme of [l|, l2|. 
The (^-function of an operator H is defined as 



(5) 



where are the eigenvalues of H and ^ represents the renormalization scale. The (^-function 
regularized effective action of the operator H is 



2ds 



(6) 



Using the proper-time method of Dowker and Critchley pj one can write the (^-function as 



Cnis) 



Jo 



1 d 
2ds 



T(s] 



(7) 
(8) 



Taking the first variation of the (^-function 

SCh 

i l^V JO 

s rco 

/ dtt'j2 

JO 



^1 dtm{6He-^-) 



'r(s] 



{n \SH\ n) 



(9) 



we obtain the regularized expectation value of the stress-energy tensor given by 



I d \ ii' 
'2ds 1 f(s) 



2 / I , 



s^O 



(10) 
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where 



T^z. [4>n{x),(f)l,{x)] 



9[X) 
2 



in -— \n) 



X] 



5H 



dgt'^ix 



(11) 



Here (f)n{x) are eigenfunctions of the operator H corresponding to the eigenvalues A„, namely, 



H(j)n = A„0„. 



(12) 



With the form of H given by ([2]) we have 



T^,u [0n(a;), 0*/(x)] = - {d^(t)*^,dy(t)n + 5i.0*,(9^0„) + gf^,. {g°'^da(t)*n,dp(t)n + 0n'n0n) 

-2c [^7^,n(0:,0„) - v^v.(0:,</)„) + (13) 



Using the Schwinger method 



2l|, Phillips and Hu j4| generalized the ^-function method 



for the consideration of the stress-energy correlators. The second variation of the (^-function 
can be written as 

s roc rco 

h^iCn = / du / dv{u + vy{uvY 

2r(s) Jo Jo 

{Tr [(5ii7)e-"^(52i/)e-^^] + Tr [{52H)e-^'' {5^H)e-^''] ] (14) 



and 



1 d \ jj,^ 



/ du dv{u + vy{uvy e' 



2ds 1 r(s) 



(15) 



Note that in this prescription an ^dditional regularization factor {uvY has been introduced. 
This is because the authors of [4| were interested in the fluctuations of the stress-energy 
tensor, that is, in the coincident limit of AT^^^,^,{x,x') where under this limit further 
divergences occur which call for an additional regularization factor. (See also [5]). However, 
our present purpose is focused on getting the correlators with two points separated, i.e., 
in the non-coincident case. Hence, apart from the fact that the expression in Eq. ( fT5l) is 
more symmetric with this factor, the keeping of this factor above is actually a matter of 
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convenience. Here we can first take the s — ?■ limit without spoiling the regularization and 
the expression in Eq. f[T5|) becomes 



1 i*co i*oo 

^ ^ n,n' 



(16) 



We shall see from the following calculations that with this expression the integrations over 
u and V effectively separate. The calculations are therefore simplified considerably. 



III. STRESS ENERGY CORRELATORS IN HOT FLAT SPACES 



We now apply this generalized (^-function method to calculate the stress energy correla- 
tors of massless, minimally coupled quantum fields in hot flat spaces. The N-dimensional 
Euclidean flat space metric has the form 



ds = dr + dx ■ dx 



(17) 



where r has periodicity /3 = 1/T. The operator H is just — □ and the corresponding 
normalized eigenfunctions are 

1 



where Un = 2nn/ [3 for n = 0, ±1, ±2, 



(2vr)(^-i)/2y^ 
. , and 



(18) 



dr I d^ ^x^l,p{T,x)(l)^^{T,x) = 5nn'S{k-k') 



(19) 



From the tensorial_ structure of the AT^j^^,^,(x, x'), one can define the scalar coefficient 



functions as follows 



23|. 



AToooo(a;,a; 
AToooi{x,x 
AToojj(x, X' 
ATQioj{x, X' 
AToijii-^x, X 



Cn 

C2l{x-X)i 

CsAj + C32(f - x')i{x - x')j 
C^lSij + C42{X - x')i{x - x')j 

C^iix - x')iSjk + C52 [5ij{x -x')k + Sik{x - x')j] 

+C53{X - x')i{x - x')j{x - x')k 



(20) 
(21) 
(22) 
(23) 

(24) 
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ATijki{x,x') = CeiSijSki + CG2{SikSji + SiiSjk) 

+Ce3 l^ijix - x')k{x - x')i + 5ki{x - x')i{x - x')j\ 
+^64 [kk{x - x)j{x - x')i + 6a{x - x )j{x - x')k 

+6jk{x - x')i{x - x')i + 6ji{x - x')i{x - x')k] 
+Ce5{x - x')i{x - x')j{x - x')k{x - x')i 



(25) 



where, owing to the homogeneity of the space, the coefficients Cab are functions of r — r' 
and |x — x'l only. 

The components of x), 0*,^,(t', x')] in hot flat space are 



^oo[0„fc(r,x),0* r,(r',f' )] 



:(^,^),C'P(^''^')] 



-UJn{(^n + (^n') -k- ik-k') 0„g(r, f ) , 0*,^, (^') ^) (26) 
-UJnk[ - Un'ki] 0„^(r, f ), 0*,^,(^'> (27) 



(28) 



-kik'j - k-kj - 6ij [uJn (w„ - UJn') + k- {k - k') 

'^nfc(^>^)'C'fc'(^''^') 

The coefficients Cab can be expressed in terms of the function /(a). 



(29) 



and its r derivatives. These expressions can be found in Appendix \M 



IV. LOW AND HIGH TEMPERATURE EXPANSIONS OF THE CORRELATORS 



To develop the low and high temperature expansions of the correlators, we look at the 
corresponding behavior of the function /(a) defined in Eq. fl29l) . In fact, for the low tem- 
perature regime, a Poisson summation formula has to be used to obtain the appropriate 
form of /(a). In these expansions we shall see that the dependence of the correlators on 
the temperature are quite different in the low and the high temperature regimes. We shall 
also discuss the corresponding singular behaviors of the correlators in the x ^ x' coincident 
limit. 
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A. High temperature expansion 



First for high temperature or small (3, ojn = 27rn//3 is always large except for n = 0. 
Hence, we need to treat the n = term separately, as follows: 



/(«) 



POD POO 

/ rfMM-"e-("^'^')'/^" + 2^cos[u;„(r-r')] / u-'e^'-e^^"^')'/^" 
Jo Jo 



4"-^r(a-l)[(f-f') ] 



n=l 



27m I X — x' I 



27™(r — r')\ / 27m\x — x' 

cos I I fsTo-i 



(30) 



Using the asymptotic expansion of the modified Bessel function Ki,{z) for large z, we can 
further expand /(a). 



/(a) = 4°-T(a- 



■/\2il-a 



1 + 



2tt 



Q-1 



r(a- 1) 



-27r|x-x'|//3 



\X — X 



/3 



cos 



An \\x — x'\J 



2ti{t - r' 



+ 



(31) 



This asymptotic series of /(«) can now be used to develop the high temperature expan- 
sions of the correlator components. They are listed in Appendix [Bl From Eqs. (IBip to 
flB6l) . we see that the density-density correlator ATq^qqq(x, x'), the density-pressure correla- 
tors ATq^q^^.(x, x') and the pressure- pressure correlators /^Tf-j^i^x^x') are of the order of 
or T^, while /^Tqqq^^x^x')^ ATg%j(x,a;') and ATq^^.;.(x, x') are suppressed by e~'^ . 

Next, we investigate the singular behaviors of the correlators in the x' ^ x coincident 
limit. Since the space is homogeneous and isotropic, we need to average over the directions 
in this limit. This amounts to replacing (x — x')i(x — x')j with 5ij{x — x'Y /{N — 1) and 
(x - x')i(x - x')j(x - x')fc(x - x')i with {6ij5ki + 6ik6ji + 6ii5jk){x - x'Y/{N - 1)^ In this 
way, the correlators are 

(Ar-l)(iV-2)r2[(iV-l)/2] 



Aroooo(a;,x') ~ 



^T^o^A^,^') 



^T^jki{^,x') 



87r^-i|x-x' 1 2^-2/32 
(iV_2)(iV-3)r2[(iV-l)/2]5, 



87r^-i|x-x' 1 2^-2^2 

(A^2 _ 7iv + 14)r2[(A^ - l)/2] 
87r^-i|x-x'|2^-2/32 



(32) 
(33) 



m-7N + 14 



(34) 
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For example, for = 4, we have 



A^^oo(^^-')-Y^(^^) (35) 



x') ~ ^ ( |f-f^|6/j2 ) (-^^i-^fc' + ^^^5,1 + (37) 

The density- density fluctuation is much larger than the other fluctuations. However, when 
N is large, all these three fluctuations have the same iV^r^[(iV — l)/2] dependence. 



B. Low temperature expansion 



To find the low temperature expansion for the correlators we again consider the function 
/(a). The definition in Eq. ( I29l) is not suitable for this purpose. An appropriate form can 
be obtained using the Poisson summation formula (See, for example, 24j|). 



n=— oo n=— oo ^ ' 

where g{k) is the Fourier transform of g{x), 

m- 



dx 



e '■^''g{x) 



Taking g{x) = e 



-Att'^Ux'^ / l3'^+i2lT{T-T')x/ 13'^ 



, we have 



E 



-UuJ^+iuJn(T-T') 



E 



'Attu 



and the function /(a) becomes 



f{a)= X: -^^[ix-xr + ir-r'-n^r] 



21 



Air 



(38) 



(39) 



(40) 



(41) 



In this form one can expand in powers of for large /3 or low temperature. 



47r 



2q-2 



(42) 
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Using this expansion we can develop the low temperature expansions of the correlators. 
They are listed in Appendix [Cl We make two observations: 1) From Eqs. (IC1|) to flC6l) . we 
see that when /3 — )■ oo or T — )■ 0, the correlators reduce to the ones in a i?^ space. One 
can check that the values we obtain here are the same as the ones in jo]. 2) The leading 
temperature dependence of the correlators are all of or . 

To investigate the singular behaviors of the correlators in the coincident limit in this low 
temperature regime, we first take the r' — t- r limit. In this limit ATq^qq^ and ATq-^, vanish. 
We then average over the directions in the other correlators giving 

N' 



^^0000 ^ ) 



N{N -l)V^{N/2) 

87r^|f-f'|27V 



4C(iV) 
A^- 1 



\x — X 



(3 



(43) 



'r\2N 



4{N^ - 3N + 4:)C{N) f\x 
(iV3 _ 4Ar2 + + 4) 



X 



N' 



^T^^o,i^,x') 



T'{N/2)S,, 



4(iV_l),.p_J-r l-4(iV-l)C(iV) 



\x — X 



/3 



N 



(44) 
(45) 



X 



8(A^- l)%^|f 
4(A^-1)(A^2^7A^-4)C(A^) f\x-x"^ ^' 



{m - 7m + 15m - a^ - 4) V /? 



{m + l)T^{N/2){6,k6,i + 6u6,k) 
4(A^- l)27r^|f -f'|2^ 



4(A^- l)C(A^) f\x-x' 



(m + 1) 



(3 



N 



(46) 

(47) 



From the signs of the temperature dependent terms we see that the finite temperature effect 
tends to increase the density-pressure fluctuation while it tends to decrease all the other 
fluctuations. 

To see their relative magnitudes we examine the four dimensional case. Setting A^ = 4 in 
the above we get 



^T^O^,i^,x') 



2tt^\x-x'\^ 



3tt*\x — X 



■/|8 



1 



1 + 



2ti^ 


^ \x 


— x'\ 


135 








( \^ 


-x'\ 


45 





(48) 
(49) 
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Here the finite temperature effect is larger for tfie current-current and tfie pressure-pressure 
fluctuations than the others. On the other hand, for large values N, the current-current 
fluctuation goes hke 4:N({N), while the other fluctuations all go like A({N)/N. 

V. CONCLUDING REMARKS 

In this paper we have calculated the stress-energy tensor correlators of a massless, mini- 
mally coupled scalar quantum field in hot N-dimensional flat spaces. With the help of the 
Poisson summation formula we are able to develop low temperature as well as high temper- 
ature expansions of these correlators. Low and high temperature regimes are determined by 
whether T is smaller or larger than \x — x'\, the only dimensional parameter in the problem. 
From these expansions we see that the correlators have rather different behaviors in these 
regimes. For low temperature the correlators all have finite temperature corrections of the 
order . However, for high temperature, the density-density, the density-pressure and the 
pressure-pressure correlators are of the order T^, while the other correlators are suppressed 
by e-^. 

We have also investigated the quantum fluctuations of the stress-energy tensor by con- 
sidering the singular behaviors of the correlators under the x' ^ x coincident limit. Here 
we first take r' ^ r and then average over the directions. In the low temperature regime we 
see that the finite temperature contributions tend to decrease the magnitude of fluctuation 
because the vacuum and the thermal fluctuations are of opposite signs, except for ATqq^j. 
On the other hand, in the high temperature regime, the unsuppressed correlators all have 
similar magnitude which go like A'"^r^[(A^ — l)/2] for large N. Note that in this limit, both 
13 = 1/T and \x — x'\ are small but \x — x'\ > (3 is always assumed in our approximation. 

It would be interesting to see how these behaviors would change with the introduction of 
another dimensional parameter. For example, we can consider massive scalars or background 
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spacetime with non-zero curvature, like in {g]. We do not expect the high temperature be- 
havior to change. However, the low temperature behavior might be altered. This makes 
careful consideration of finite yet small temperature corrections to established results in- 
volving quantum vacuum quantities, such as negative energy density and theorems derived 
therefrom, such as the quantum interest principle, worth its while. We hope to address 
these foundational issues of importance and report our findings certainly on, if not before, 
the occasion of Professor Dowker's 80th birthday celebration. 
Here's to wishing you, Stuart, many happy returns! 
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Appendix A: CoefRcient functions 



The coefficient functions of the stress-energy correlators in Eqs. (120|) to (125|1 : 



Cii(t — t', \x — x' 

1 



2(47r 



+ (iV-l) 



1 



[x — x'Y 



N -I 



(x — x')^ 



Vh-2)./,^.1,/(^H-3) 



C2i(7' - r', |x - x'l) 
1 



2(47r)^-i/32 



X 



2' 
iV + 1 



A^-l 



+ 



[X — X' 



iV - 1 
f{^— + 2)drf{ 



2 

A^- 1 



1) 



1 iV-1 1 N-l N-l 

+ 2/(^- + + 2) + 2^(^- + 3)5./(^-; 
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C3l(T-TMf-f|) 
1 



X 



2(47r)^-i/32 
"iV2 -iV-2 „.7V- 1 



-/(- 



'l\2 



{x — x') 



N -I 



+ + 1) - - + mn-. 



.N-l. 



N-1 



2 +2)a^/(^) + 2/( ^ 



N-l 



+ 1) 



TV - 1 iV - 1 

-(2Ar+l)/(^— + 2)/(^— + 1) 



[x — x') 



•/\4 



2 

iV- 1 



+ 3)/(^ + l, + /(^^ + 2)/(i^ + 2) 



C32(t - T , |f - f I) 
1 



2(47r) 



Ar-l/92 



X 



(x — x') 



•/\2 



fi 



N-l 



+ 3)/( 



N-l 



+ !) + /( 



N-l 



+ 2)/( 



iV- 1 



+ 2) 



C4i(r — t' . \x — x'\) 



C'42('?' — 7"', |a; — if 
1 



2(47r)^-i/3 



X 



1 TV 
- + W(- 
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C5l(T-TMf-f|) 
1 



2(47r)^-i/32 



X 



i±i/(^+iR/(^+i) 



,iV- 1, 

y 



•/\2 



{x — x') 



/( 



iV- 1 



2 

iV-1 iV-1 
+ 2)9./(^- + 1) + 2/(^- + 3)9./(^-) 



1 - 1 N -I 



C52(T-r', |f-f'|) 

1 



2(47r)^-i;g2 1^2-' ^ 2 



C53(r-r', |f-f'|) 
1 



2(47r)^-i/52 1^ 2-' ' 2 



1 iV - 1 iV - 1 

-o/(^^ + 2R/(^^ + l) 



(A7) 
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Appendix B: High temperature expansion 



Using Eq. ( 13T|) we have the high temperature expansions of the correlators 
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Appendix C: Low temperature expansion 



Using Eq. ( H2|) we have the low temperature expansions of the correlators 
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